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INTRODUCTION 


This document contains r eference materials, calculations, and trade studies used in the 
analysis and selection of YSTAR components. The document is organized by each major YSTAR 
system each of which contains material pertinent to that system. Many sections contain simple 
graphs and tables used to make qualitative comparisons of various YSTAR component candidates. 
Equations and/or calculations used for a particular analysis are also Included where applicable. 
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Mission Plans 
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aY analysis 


SLR Mission 

This analysis determines the Impulse maneuvers required for executing LEO to 6E0 
Hohmann transfers for the SLR mission. Two possible methods of transfer are considered for 
this orbit change. The first aV calculation applies to the LEO-GEO transfer involving the 
transfer from LEO to a high-earth orbit followed by a plane change to place the vehicle in 
geosynchronous orbit The total aV required for the method 1 transfer is found to be 
approximately 18,040 ft/sec. The sec&w calculation is fGr a LEO-GEO transfer in which the 
orbital plane change Is executed by a two-impulse maneuver within the altitude change. The 
total aV required for the method 2 transfer is found to be approximately 14,271 ft/sec. 

■ A 

IYA Mission 

This analysis determines the impulse maneuvers required for executing the 53.27* 
plane change in low-earth orbit (28.5* Space Station inclination plus 23.27* Earth inclination 
relative to the eclliptfc plane). The total aV required for the constant orbit plane change is 
approximately 22,504 ft/sec. 

The calculations associated with this analysis assume a low-earth orbit altitude of 230 
miles and a high-eerth orbit altitude of 22,000 miles. Derivation of equations used in aV 
calculations is given in (Ref. 4 1 ). 
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Table i - Electric Propulsion System Comparisons («**• 


Ada L 


AdK high tap *mp*dty 


» W ■ 

smpaacy 
highly developed 


Nghtap 
high thrust 


Otecomplex electrodesrosion t*mp.Bmit*d power efficiency highcurrent 


high powsr 


Table 2 - Liquid Chemical Rocket Engine Comparisons (f *6 is) 


--A ^-—-< n I *1 I ^- A - 1 < - i - ! gJ - ^-•- 

nocNMoyniMmWUnn innMfvviGWcn^ni. 


lap-482 secs 

Engine length (retracted)-) .02 m 

Props*** LH2.LOX 

(extanded)-3.91 m 

Thrust ma»-68750 N 

Node ext dfemetarel .96 m 

min-2225 N 

Area redo-1300:1 

MR(OF)-aort 

Engine mass-196 kg 

' • p . 

—-* |f4 |- -* ^ —--— % -i _ ^_ i __ 

nw ot rmmwf wwinciq uwa iranpirvinKM engm: 

ap—w secs 

Engine length (retrected)-l.02 m 

Proped* LH2.LQX 

(«xtended)-3.05 m 

' Thrust mas^66750 N 

Nod* exftdtometer-1.63 m 

min-2225N 

Area rado-640:l 

MR(OP)-6.0rt 

Engine mass-204 Kg 


Aerojet Afranced Orbital Transfer 
lap-482 SSCS 
P rops- * LH2.LOX 
Thrust mix* 13850 N 
m ir w44S N 
MR(OF)-QlO:1 


Engine: 

Engine Isngth (retract*d)-0.98 m 
(extende d )-1.96 m 
Nod* ext d a mfr-0.76 m 
Are* ratio-1200:1 
Engin* mass-57 Kg 


Prat and Whitney RL*10 Category IV: 
lap-470 sacs 
PropatfanMJg.LQX 
Thrust max-68750 N 
min-16688 N 
MR (OP)-6.0:1 


si 


Engin* length (retracted)-1.45 m 
(*xtend*d)-2.9G m 


Nods exit d ia meterel.68 m 
Area ratio-401:125 
Engin* mass-193 Kg 
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Table 3 - Dual Fuel / Milted Mode Rocket Engine Characteristics <o 


Pa raw ter 


Mode 1 

(LOX/RP-l/LHj) 


Both 

Modes 


Mode 2 

(lox/lh 2 ) 


Vacuua Thrust (Ibf) 20.000. 

Vacuus Delivered 

Specific Impulse (sec.) 418.6 

Chamber Pressure (psla) 2000. 

Total Propellant Flowrate 

(1b/sec) 47.8 

Mixture Ratio 4.25 

Nozzle Area Ratio 

Nozzle Exit Diameter (In.) 

Engine Length (in.) 

Total Engine Weight (1b.) 


9915. 

460.6 

1007. 

21.5 

7.0 

400. 

49.3 

t'zm.z* 

557. 


° 2‘"2 


_ SINGLE-WOE Oj/!H OTV'S 

MIXTURE HAT 10 -6.1 ~ 


Oj-COOLED ENGINE; 



DEPLOY P/L - 7926 SAM LB 

LB 


9476 LB 



DRY NT. - 5150 LB 


DRY NT. - 5061 LB 


DRY Iff. - 4744 LB 




ItlXED-MOOE o 2 > , wk/h 2 OTV 

MX. DEPLOY PAYLOAD 

DEPLOY P/L • 9601 IB 

* 



DRY Iff. • 4469 LB 


DUAL-FUEL P!3. 

NON. DEPLOY PAYLOAD 
7926 LB 



DRY Iff. • 3925 LB 


Oj/fflH OTV 

O 2 -COOLED ENGINE 
6093 LB 



DRY NT. • 3491 LB 


OTgARJSON.QFSINGLE-MODE ANDNl 


IGNITION NT. WITH MX. PA 


SBnUflfti 


Figure I - Single Us. Mined Mode OTU Comparisons (Ref. b) 
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Thermal Protection 

The selection of the type of radiator was based on survivability, and the membrane area 
per Btu/h output Typically the higher the area required the higher the weight of the device. The 
estimated survivability was based on research. For VSTAR the heat pipes are to be used as thermal 
buses between fluid systems and the choice of these based on the same characteristics as the 
radiator. The data used to determine the ft*2/Btu/h comes from several sources. The Btu/h is 
converted from Watts and this wattage Is the maximum output of the device under choked flow 
conditions. As can be seen in Table , the best heat pipe for the thermal bus is the axial ly grooved 
variable-conductance heat pipe, but as a radiator, Its survivability Is far to low. As for the 
radiator, the rotating bubble membrane radiator is the superior choice as n is highly survlvable 
and does not waste liquid. (Ref. °\, 

» 

Table W. Radiator and Heat Pump Data 


Jm 

ft*2/Btu/h 

Survivability 

Osmotic Heat Pipe 

0.00293 

• 

medium-low 

Heat Pipe Sandwich Panel 

0.00007 

medium 

Double Walled Artery Heat Pipe 

3.264E-7 

medium-low 

Axially Grooved Variable- 



Conductance Heat Pipe 

1.3I0E-8 

medium 

RBMR 

0.00019 

high 

Liquid Droplet Radiator 

NA* 

medium-high 

Heat Pipe as Radiator (ave) 

0.00075 

medium to low 


* this device sprays the hot liquid directly into space but would weigh the same as the RMBR 
before the liquid was used up. 


12. 















Material Selection 

Table Sshows the stress allowable and density of each material that was considered. The 
exact cost of each material produced is not known, however, based on research estimated relative 
costs have been determined. Compatibility Is known and Is accounted for in the design. The 
material with the smallest densfty/stress is the strongest per unit mass. Metals withstand debris 
and micrometeoroid Impact much better than do composites with epoxy bases, because metals tend 
to buckle before they are punctured and epoxys tend to shatter. 

Therefore a metal is preferred for debris protection. Careful consideration of this data with the 
micrometeoroid protection results in the choices made for the vehicle components.(Ref. n& 29 
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Table 5. Material Properties (Ref. £9 


material 


density. 


SSL 


(psl) 

(1bm/1n A 3) 



4130 Steel 

70000 

0.283 

4.0428E-6 

low 

2024 Aluminum 

42000 

0.100 

2.3809E-6 

very low 

Boron Epoxy 

85333 

0.072 

8.496 IE-7 

medium 

Graphite Epoxy 

100000 

0.057 

5.7000E-7 

medium 

Aluminum Oxide/ 





Aluminum 

42222 

0.100 

2.3684E-6 

low 

Boron Aluminum 





2024 

100000 

0.100 

1.00E-6 

high 

* 

Boron Aluminum 





6061 

96889 

0.100 

1.032E-6 

high 

Graphite Aluminum 





GY70 2024 

37911 

0.100 

2.638E-6 

high 

Graphite Aluminum 





V60054 2024 

43156 

0.100 

2.317E-6 

high 

Graphite Aluminum 





T50 2024 

31111 

0.100 

3.214E-6 

high 


1+ 




















Debris and Micrometeoroid Protection 

BUMPER is a program provided to us by NASA/Marshal to determine the chance of 
penetration of a vessel BUMPER uses the equations based on some simplifying assumptions and 
your inputs. There are two modes corresponding to single wall and double wall vessels. For double 
wall, inputs are surface area, sheild thickness, sheild stand-off distance, vessel thickness, and 
whether or not MU is being used. For the single wall surface area, vessel thickness and MU use is 
inputted. The program then returns the percented chance of debris or micrometeoroid protection, 

(Ref. u). 
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Table £ EC/LSS System Functions 10) 


ecus FUNCTION 

. MAJOH fOUIAMCNT 

• ATMOOAMCNC HNCSSUNC A COMHOSITION CONTROL 

• MNC OCTCCTION 1SUAARCSSION 

MlSUil RCGULATION 

FOdTAMI OXVCSN S VST CM 

smok«/si*« octsctons 

TIN* SURRAS SSI ON SVSTtM 

• MODULI T|MMI9ATV«I 1 HUMIDITY CONTROL 

OCHUMIOtfICAtlON VfNTILATiON FANS 

AIN COOLING MCAT t XCHAMGC NS 

• ATMOS***** *CVtTALlSATlOM 

• CO) eONTNOL/OSMOVAL/NCOUCTION 

• 0 |*N|MAK|UP 

• TNACS GAS MOMtTOAlMO 1 CONTROL 

CARSON OIORIOC RSMOVAL ANO ORYOCN OCNSRATION 

eOLLICTION (MCROCNCY OXT 6 CN ANO 

CAHSON OIORIOC REDUCTION NlTROCSN JTORAOC 

contamination control 

OOON CONTROL 

ATMOSAMSRS MONITOR MO 

• WATCH MANA4CMCNT 

• NASTI WATCH COLLCCTlQNWNOCMBiNO 

• WATCH QUALITY MONITOR MS 

• trOHATC A OfSTRIOUTION OA HCCOVCNCO WATCH 

SVAAONATtON TUNIS MCA TIOM 

MATIN DUALITY MONITONING 

MATIN ST ON AGS 

• NASTI MANAGS Mf NT 

• C0LLSC***O£SSS UNfN* 

• COLLfCUSTONt MCAL MAlTtfl 

MASTS eOLLICTION ANO STONAGC 

CMSNGfNCV WAITS COLLICTION 

HOT/COLO WATS* SUTRLV 

• CVA SUAAOHT 

• AHOVIOC CSACNOAALCS/SUAAQHT TO CMU A MNW 

• AHOVIOC LI AC SUAAOHT SSRYICU TO 

aihlocx/mvacaaahic hacility 

»* 

• 

sum ANO SAexTACXS 

NfGNANGf STATIONS 

AIN LOCK SUNONT 



* 




Table 1 • EC/LSS Performance Requirements (jut. \o) 


parameter 

UNITS 

operational - 

OCCRAOEO m 

CO* NAWTIAL MISS 
TSMAfNATUNS 

MMHQ 

AO MAX 

7.1 MAX 

0C6I 

•*7S 

10*41 

OSW TOIMT 12) 

OfO H 

40*10 

31*70 

FOTAOLS WATS* 

u/mam-oav 

IS-S.1 

AS (31 

MVCtSNS WATS* 

LS/MAN—OAY 

IS (St 

1(31 

WASH WATER 

ls/man-oay 

71(31 

14(31 

VINTILATION 

AT/MM 

11*40 

10*100 

Of AAHTIAL 
>REMUR« (41 

nuk 

A7-A2 

2.4-34 

TOTAL HNCSSUNC (SI 

NU 

107 O* 14.7 

10-2 0* 14.7 

OILUTt OAS _ 

TRACS CONTAMINANT (Si 

MOM*. 

3o 

$o 

MIC *0*0 NO AN ISMS 

ctu/hH (It 

SOD 171 

710 <71 


WOTS* 

111 

at 


ut 


lit 

lit 

171 

<•> 




OKONAOSO LfVCLS M«*T “TAIL OAfNATtOMAi.** CRITERIA. 

mutivi MUMioiry shall m within thc ranoi oh 3-75 hscint. 

MINIMUM. - 

IN NO CAM SMALL TM| O 7 AANTIAL MUSUM M MLOW 2.3 ASIA. ON 
THC Oj CONCI NT NATION CRCCIO SJ AC RCINT OH THC TOTAL ARCSSUNC AT 
M.7 h4a ON 30 AtHCXNT OH THC TOTAL HNCSSUNC AT 10 J. 

ALL SYSTIMS SMALL M COMAATISLC WITH BOTH 10 J A NO 14.7 ASIA TOTAL 

ANIMUNC. 

CHU - COLONY HONMINO UNITS. 

TMCSS VALUCS NCHLCCT A LIMITED SAM. MO WIOCLY SANCTIONCO STANOAHOS 
ANC AVAILASLS. 

■AMO ON NH« RO40.1S. I. 
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Table % - FC/LSS Ruerage Design Loads •«& 


MSTAOOUCO* 
LSAKAOS AM 
IVAO} • 
•VACO* 
MSTAOODC CO} 


PAAATKSN HjO 

mooOcti 




• oammch*o 

- FOOO Rfttr 

- MSTASOLIC HjO 

• CLOTHS WASH HyO 

• HAMO WASH HiO 

- SH O W S H^O 

- SVA HjO 

- OSR3P1RATION AMO ftSSWRATION 

- URINAL SLUSH MjO 

- UHMM HjO 

- sooosoiu 

- A 000 H 70 
OACI 


- UHIMS SOLIOS 

- MCAL SOLIOS 

- SWSATSOUOS 

- SVA WASTSWATSH 

- CHAR CO AL WSOUISSO 

- MSTASOLIC SSNSISLS MSAT 

- HVOISMS LATSMT HjO 

- ROQO RRSPARATtON LATSMT HjO 

- LAUNORV LATSMT H*0 

- WASH HjO SOLIOS 

- SHOWK/MAMO WASH H*0 SOLIOS 

- AIM LOCK OAS LOSS 

- TKASH 


h*o 


- THASH volums 


1A4 LS/MAM OAV 
LM LS/DAV TOTAL 
1.22 LS/S MS SVA 
IAS LS/S HA SVA 
US LS/MAM OAV 
4JM LSAAAN OAV 
1JSLS/MAN OAV 
0.7S LS/MAM OAV 
27JS LS/MAM OAV 
4JM LS/MAM OAV 
MS LS/MAM OAV 
%M LS/S HM SVA 
4JB LS/MAM OAV 
IAS LS/MAM OAV 
&31 LS/MAM OAV 
IAS LS/MAM OAV 
IAS LS/MAM OAV 
IAS LS/MAM OAV 
#.13 IA/MAM OAV 
047 LS/MAM OAV 
SA# LS/MAM OAV 
2AD LS/S HU SVA 
0.13 LS/MAM OAV 
7SOS STVI/MAN OAV 
OAS LS/MAM OAV 
OAS LS/MAM OAV 
0.13 LS/MAM OAV 
0.44ft 
0.12ft 

1A3LSS/USS 
1A0 LS/MAM OAV 
0.M 0T4/MAN OAV 


Table R - EC/LSS Technology Requirements Ou*. 10) 


- SSCAL MASTS MANAOSMSNT 

- TRASH/OOOO MAMAOSMSMT 

- SSMSOM OS VS LOOMSMT 

O MASS OAUOIMO 


•AJR/WATSA QUALITY 

- WAT8R RSCLAMAnON/OflOCSSSINO 3VSTSMS 

- RSOSMSRATIVS CO] RSMOVAL/MSOUCTIOM SV3TSM 

- MARS ATMO S DIISRS OROCfSSNM SVSTSM SOM OXVQSM. 
MITROOSN A WATSR 
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FAMILY OF PRIMARY BATTERY CELLS 
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COST ANALYSIS 


21 





















There are several equations required in cost analysis. These equations are given below for 
each of the subsections, vehicle costs and program management costs. Vehicle costs are in 1978 dollars 
which must be multipled by 2 to approximate f 988 dollars. The management costs are functions of the 
vehicle costs. C***. iO 

Vehicle Costs 

DDT&E 

Structural Cost * (2)( 1.8531)(Structural Weight) 0 -491 

Electrical System Cost - (2X0.6000X Electrical System Weight) 0,584 

Propulsion System Cost * (2X0.1075X Propulsion System Weight) 0,878 

» 

* 

FH 

Structural Cost * (2)(0.4445XStructural Weight) 0 - 440 

Electrical System Cost«(2)(0.0422XElectrical System Weight) 0,784 

Propulsion System Cost»(2X0.1158)(Propulsion System Weight) 0,550 

Management Costs 

Integration, Assembly & Checkout 
DDT&E-9.1208(STE)°- 4 °1 
FH - (Total FH) 0,832 

System Test and Evaluation 

DDT&E - 13.101 (Total FH) 0,672 

FH is zero. It is accounted for in cost of operation. 






















S/stems Engineering & Integration 
DDT&E « 0.4519(DI5G) 0 - 876 
FH « 0.3750< Total FH) 0 - 855 

Program Management 
DDT&E - 0.6952( DIS6S) 0 * 731 
FH - 0.3146(FIS) 0 * 798 

Suplementary Equations 

STE - 11.232(Total FH) 0 - 672 

DIS0 * Total DDT&E + Integ., Assem. & Check (DDT&E) + Sys. Test & Eval. (DDT&E) 

DISGS - STE + Total DDT&E ♦ Integ., Assem. & Check (DDT&E) ♦ Sys. Eng. & Integ. (DDT&E) 
FIS » Total FH ♦ Integ., Assem. & Check (FH) + Sys. Eng. & Eval. (FH) 
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Optimization Methods 


The following describes the derivation and assumptions made for each of the mathematical 
expressions describing the calculated parameters. The analysis is based on a division of the 
configuration Into the components shown in figure 3. 

The mass of components In many cases ars tn& -jm of the masses of the subcomponents. 

Their derivations are discussed under the appropriate sections and are not repeated here. 
Calculation of m^n: The mass of the cabin is constant and not dependent on the 

acceleration loads during flight This is justified by assuming that It will be launched in one piece 
from the space shuttle and will necessarily encounter higher accelerations during launch than 
those experienced for a typical mission. 

m cab1n * m crew + m air + m $tructure + m shie1d1ng + mother 
-600+ 150+ 1140+ 157 + 200 
- 2247 


Calculation for mw e support : mass of the life support system includes the mass of 

the power supply system and is equal to a base mass plus the mass of the consumables (breathing 
oxygen, water, etc) used per hour by the crew. 


mitfe support * f (* of crew, Tg^) 

» 1075+ 1.39 nT^ 

Calculation of mgargo module : Tt1e mass of the cargo module is dependent upon the mass of 
structures above it which it will have to accelerate during thrusting periods, as shown in figure 5 
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Figure 13- Idealized Components for Rnolysis Purposes 




















m cargo module * f ( m cab1n» m 11fe support* m pay> ®max) 


An approximate function for the mass can be developed by considering figure t f. 

Let A be the cross sectional area of a column of density, r, and maximum allowable stress, 
o, and length 1. The material must withstand the stress generated from mass m and acceleration a. 
The force exerted on the column is: 

F«ma 


and Is dqual: 


F-Ao 

« 

Equating and solving for A gives: 

A«ma/o 

The mass of the column is then: 

m column * A1 r 

*r1ma/o . (l) 

This equation is based on the assumption that the mass of the column does not contribute 
any to the stresses in It. This is a reasonable assumption if m » mcolumn. 

An extrapolation can be made from the column of Figure 6 to the proposed truss structure 

of the cargo module in Figure s if the additional assumptions are made: 
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]) The cross sectional area Is equal to the sum of the cross sectional areas of the 
constituent members of the structure. 

2) The maximum allowable stress In the supporting members Is much lower than the 
compression yield stress to account for lower buckling stresses In a slender column. 

3) The effective density is higher than the actual density to account for the presence of 
cross beams and diagonal members which do not bear the compressive loads. 

Typical r/o values for steel and aluminum are 3.8 x 10"^ and 2.9 x 10“® [ibm In /Ibf], 
respectively. For Aluminum with the maximum allowable stress to be 1 /10 the compressive 
/ield stress, the effective density to be 3 times the actual density, and a length of 20 feet, the 
following expression results: 


^support structure * 0*® 1 m 9max 


Including an approximate value for the mass of the metlorlte shielding gives: 
msupport structure *150 + 0.017 m amax 
Where: 


m * m cabin + m hfe support + m pay 

Calculation of m^ : Two possible configurations for the arrangement of the fuel tanks 
are shown in figure 5 . 

The talks in Figures*. are of equal radii, but because of different volume requirements 
they are of different lengths. They employ a monocoque structure which resists internal pressure 
and the forces due to acceleration and the mass above it. Figure ?«shows two spherical tanks with 
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a separate support structue for resisting the acceleration loads. The tanks in Figure are 

efficient in that they provide a dual role In resisting the Internal pressure and the loads due to 

• » 

acceleration. However, the nonsherical geometry is Inefficient in providing a light tank for a 

given volume. Additionally, the monocoque skirt between the two tanks will provide an excellent 
path for heat conduction between the two tanks at different temperatures, which if carrying LO 2 

and LH 2 will be over 100 degrees ranklne difference. This will complicate thermal control of the 

two tanks and Increase bolloff losses. The arrangement of Figure78 was chosen because it allows 
the use of spherical tanks which weigh less for a given tank volume, and since they are not actively 
involved in resisting the acceleration loads, they can be more completley thermally isolated. 

The mass of the tank can be expressed In the following form: 

m tank * f ( m fuel* r fuel» m evap» Pfoel* ^min* r metal) 


The required volume is given by: 

v tank * ( m fuei + m evap ) / r fuel 
The radius is given by: 

rtank-f 3 v^/ 4 pi) 1/3 


From fluid statics the pressure of a liquid of density, r, under an acceleration, a, at a 
depth of h Is: 

P * P 0 + r a h 
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Applying this to a spherical fuel tank.: 


Pbottom * Pvapor + 2 rf W i A^p 

It can be shown that for a constant thrust, the greatest pressure will occur at the bottom 
of the fully filled fuel tank even though the acceleration Is at a minimum. Assuming a thin walled 
pressure vessel the thickness required Is: 

ttank * Pbottom r tank / 2 o 

* 

The area Is: 

Afenk * 4 r tank ^ 

The volume of metal In the tank Is: 


v metal In tank * *tank ttank 

And the mass of the tank Is: 

m tank * y metal In tank r metal 

Calculation of m^op: The fuels are cryogenic and are assumed to be continual !y receiving 
a heat input which holds them at the boiling point for a given pressure. Any heat input is used in 
vaporizing the liquid into its gaseous state, which is immediatiy vented, thus keeping the pressure 
constant. 
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The particular Insulation chosen conducts heat thirty times faster In the direction parallel 
to the layers than In the direction perpendicular to the layers. It will be assumed that the heat 
received from radiation on the Illuminated side of the tank will be conducted to the unlllumlnated 
side allowing an average surface temperature Tout to 1)6 used in the below equations. 

It Is assumed that this is the only source of heat input Into the tanks for the entire 
mission. This Is not true however since there will be a radiated and conducted heat transfer during 
the firing of the engines due to the fuel tanks being in close proximity to the hot engine and exhaust 
plume. These calculations may be dependent on many variables, but only the heat due to the 
radiation on the Illuminated side of the tank Is considered In the following calculation. 

The heat leak rate Into the tank Is given by: 

Qdot" k tn$ul A tank (tout “tin V tinsul 

Then the total heat leek is simply: 

Qtotal * Qdot ^end 

And the mass evaporated is: 

m evap * Qtotal 1 ttyapor 

Calculation of min$ui: The insulation is assumed constant thickness and covers the entire 
spherical tank. 


m insul * r 1n$ul A tank tinsul 
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Calculation of msupport structure The equation for the mass of the support structure Is 
developed along the same lines as m^go module R Is assumed to be of the form: 


^support structure * 200 ♦ 0.017 m amax 


Where: 


m * m cab1n + m 11fe support + m pay + m tank 
200 - approximate mass of meterolte shielding 

Calculation of mgngine : The mass of the engines are assumed to be known for a given 

configuration and are taken from the available literature. 

Calculation of mb rnou t 80(1 m total : The burnout mass Is simply the sum of all the masses 

calculated above. The total mass is the sum of the burnout mass and the mass of the fuel. 

Calculation of amax and amin : The thrust Is assumed constant over the burn time. By 

assuming specific thrust vs. time profiles, the maximum accelerations could be reduced while still 
allowing relatively short burn times, however, here It Is assumed they are constant. 

The accelerations then are simply: 

amax * Thrust / nntjrnout 


amin * Thrust / m^ai 


Calculation of mass ratio: It Is assumed their will be no losses due to gravitational and 
drag effects. The mass ratio required then for a given aV and ISP Is: 
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R » Q A Y / Q ISP 


Calculation of mf^i: 

From the definition of mass ratio: 

R * ( m^rnout + mfuei) / m brnout 
Solving for m^i gives: 

mfuei * fflfcpnout ^ “ 1 ) 

There is one other assumption made in the program that needs to be pointed out. The boiloff occurs 
before the engines ignite. This results in a higher mass ratio than actually needed since in reality 
the boiloff would occur continually and some of the fuel mass would not be imparted the full 
29,000 ft/sec av. 

Calculation of burn times: 

For a given aV: 

mfinal * mental e’ AV/ 9 |SP 

The mass of the fuel used Is: 

mfuei * minitial - mfinal 
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The thrust can be expressed as: 

Thrust * mfot fuel vexhaust 
Since the values In the above equation are constant 
Thrust ■ mfuel g ISP/1 
Solving for t 

t« mfuel g ISP / Thrust 
Correctina for units: 

t * mfuel g ISP / Thrust ge 
t ■ mfuel ISP / Thrust 

The program listing and program runs are included below. Following this are graphs 
showing data on the liquid hydrogen and liguld oxygen as used in the program. 
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Ofn*M*iAiT\otA 9«o feQAM Output 


0XY6EN/HYDR0GEN RATIO t. 3 : 

7 



OXYGEN/RPJ RATIO 

E. 3 t 

4. 23 



ISP 

E. 3 * 

420 



ACCELERATION - MAX. 

Eg'a3 t 

2. 42 



ACCELERATION - MIN. 

Eg'a3 t 

0. 28 



DELTA V - OBTAINABLE 

Eft/sac3 t 

29000 



MASS RATIO 

C. 3 t 

8.336368 



STUCTUAL COEFF. 

C. 3 : 

0. 0733 



PAYLOAD COEFFICIENT 

C. 3 t 

0. 0496 



TIME, ENDURANCE 

Chours3 i 

72 



THRUST 

Elbf 3 t 

30000 



HASSES OF COMPONENTS 

t 




CABIN 

Elba 3 t 

2247 



LIFE SUPPORT HODUL& 

Elba 3 : 

1376.86 



CARGO MODULE 

Elba 3 t 

304.8966 



CARGO 

Elba 3 t 

3000 



SUPPORT STRUCTURE 

Elba3 : 

642.4102 



HYDROGEN TANK 

Elba 3 i 

936.6639 



OXYGEN TANK 

4 

Elba 3 s 

448.6059 



RPJ TANK 

Elba] t 

74.83677 



ENGINE 

Elba 3 t 

1000 



TOTAL BURNOUT MASS 

Elba] : 

12392. 82 



MASS OF FUEL 

Elba] t 

93396. 84 



TOTAL MASS 

Elba] t 

103789. 7 





HYDROGEN 

OXYGEN 

RPJ 

TANK - MASS 

Elba] i 

936.66 

448.61 

74.84 

- MAX. STRESS 

Elbf/in A 23: 

33000.00 

33000.00 

33000. 00 

- METAL DNSTY 

Elba/£t A 33: 

489.00 

489.00 

489.00 

- PRESS. , BOT. 

Elb£/in A 23: 

20. 14 

21. 67 

21. 28 

- RADIUS 

Eft 3 t 

8.09 

6.18 

3. 42 

- SURFACE AREA 

Eft A 23 t 

822. 78 

479.65 

147.12 

- THICKNESS 

Ein3 t 

0.0279 

0.0230 

0. 0125 

- VOLUME 

Eft A 33 : 

2219.22 

987.80 

167.79 

- VOLUME, MTL 

Eft A 33 : 

1.92 

0. 92 

0. 15 

PROP.- VAPOR PRS 

Elbf/ln A 23: 

20. 00 

20.00 

20.00 

- DENSITY 

Elba/ft A 33: 

4.38 

68.64 

93. 11 

- MASS 

Elba] { 

9681.38 

67769.66 

13945.80 

- TEMP. 

E rankina3 t 

38.38 

166.41 

58. 33 

- HEAT.VAPR2 

Ebtu/lba] t 

188.16 

91.62 

100.00 

SURFACE TEMPERATUE 

Erankln*] : 

600.00 

600.00 

600.00 

INSULATION - RHO Elba/ft A 2-in3 : 

0.20 

0.20 

0. 20 

- HASS 

Elba] : 

41. 14 

23. 98 

7.36 

- THCKNS 

Ein] : 

0. 23 

0. 23 

0. 25 

-K Ebtu-in/hr-£t A 21 s 

0.000036 

0.000036 

0.000056 

HEAT LEAK * RATE 

Ebtu/hour3: 

103.51 

46.59 

17. 85 

- TOT 

E btu 3 : 

7452.54 

3354.22 

1285.23 

MASS EVAPORATED 

Elba] : 

39. 61 

36.61 

12. 85 
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OXYGEN/HYDROGEN RATIO C. 3 

: 7 


ISP 

C. 3 

: 460 


ACCELERATION - MAX. 

Cg'a3 

t 2.45 


ACCELERATION * MIN. 

Cg's3 

: 0.35 

a ■ ' 

DELTA V - OBTAINABLE 

Cft/S»c3 

: 29000 


MASS RATIO 

C. 3 

: 7.084234 


STUCTUAL COEFF. 

C. 3 

: 0.0887 


PAYLOAD COEFFICIENT 

C. 3 

: 0.0611 


TIME, ENDURANCE 

C hours 3 

J 72 


THRUST 

Clbf 3 

: 30000 


MASSES OF COMPONENTS 

• 

• 



CABIN 

Clbm 3 

2 2247 


LIFE SUPPORT MODULE 

Clbm 3 

: 1376.86 


CARGO MODULE 

Clbra3 

2 508.9224 


CARGO 

Elba] 

2 5000 


SUPPORT STRUCTURE 

Clbm 3 

: 641.6741 


HYDROGEN TANK 

Clbm 3 

: 902.9588 


OXYGEN TANK 

Clbm 3 

t 438.6012 


ENGINE 

Clbm 3 

: 1000 


TOTAL BURNOUT MASS 

Clbm 3 

: 12253.81 


MASS OF FUEL 

Clbm 3 

: 74555.06 


TOTAL MASS 

* 

Clbm 3 

2 86808.88 




HYDROGEN 

OXYGEN 

TANK - MASS 

Clbm 3 

902.96 

438.60 

• MAX. STRESS 

Clbl/ln A 23 

35000.00 

35000.00 

- METAL DNSTY 

Clbm/ft. A 33 

489.00 

489.00 

• PRESS. ,BOT. 

Clbf/in A 23 

20. 17 

22. 01 

- RADIUS 

Cft3 

7.99 

6. 10 

• SURFACE AREA 

Cft A 23 

802.17 

467.62 

- THICKNESS 

[in 3 

0. 0276 

0.0230 

- VOLUME 

Cft A 33 

2136.36 

950.87 

- VOLUME, MTL 

Cft A 33 

1.85 

0. 90 

PROP.- VAPOR PRS 

Clbf/in A 23 

20. 00 

20. OO 

- DENSITY 

Clbm/ft A 33s 4.38 

68. 64 

- MASS 

Clbm 3 

9319.38 

65235.68 

- TEMP. 

Crankinm3 

38.38 

166.41 

- HEAT.VAPRZ 

Cbtu/lbm3 

188.16 

91. 62 

SURFACE TEMPERATUE 

Cranklnm3 

600.00 

600.00 

INSULATION - RHO Clbm/ft A 2-in3 ! 

0. 20 

0. 20 

- MASS 

Clbm3 : 

40. 11 

23. 38 

- THCKNS 

tin] j 

0. 25 

0. 25 

-K Cbtu-in/hr-it A 23i 

0.000056 

0.000056 

HEAT LEAK - RATE 

Cbtu/hour3 

100.91 

45. 42 

- TOT 

Cbtu3 

7265.84 

3270.10 

MASS EVAPORATED 

Clbm 3 

38. 62 

35. 69 
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OXYGEN/HYDROGEN RATIO C. 3 


ISP 

C. 3 

ACCELERATION - MAX. 

Cg's] 

ACCELERATION - MIN. 

Cg'«3 

DELTA V - OBTAINABLE 

Cft/s*c3 

MASS RATIO 

C. 3 

STUCTUAL COEFF. 

C. 3 

PAYLOAD COEFFICIENT 

C. 3 

TIME, ENDURANCE 

Chours3 

THRUST 

Clbf 3 

MASSES OF COMPONENTS 

« 

« 

CABIN 

Clbm3 

LIFE SUPPORT MODULE 

Clbm] 

CARGO MODULE 

ClbmJ 

CARGO 

Clbm 3 

SUPPORT STRUCTURE 

Clbm 3 

HYDROGEN TANK 

Clbm3 

OXYGEN TANK 

Clbm 3 

RPJ TANK 

Clbm 3 

ENGINE 

Clbm 3 

TOTAL BURNOUT MASS 

Clbm 3 

MASS OF FUEL 

Clbm] 

TOTAL MASS 

Clbm3 


TANK 

- 

MASS 

Clbm] 


- 

MAX. STRESS 

Clbf/in A 23 


- 

METAL DNSTY 

Clbm/ft A 33 


- 

PRESS. , BOT. 

Clbf/in A 23 


- 

RADIUS 

Cft3 


- 

SURFACE AREA 

Cft A 23 


- 

THICKNESS 

C in] 


- 

VOLUME 

£ft A 33 


- 

VOLUME, MTL 

Cft A 33 

PROP. 

- 

VAPOR PRS 

Clbf/in A 23 


- 

DENSITY 

Clbm/ft A 33 


- 

MASS 

Clbm] 


- 

TEMP. 

Crankin«3 


- 

HEAT.VAPR2 

Cbtu/lbm] 

SURFACE 

TEMPERATUE 

[rankine] 

INSULATION - RHO Clbm/ft A 2-in3 



- MASS 

C lbm] 



- THCKNS 

Cin 3 



-K Cbtu-in/hr-ft A 23 

HEAT 

LEAK - RATE 

Cbtu/hour] 



- TOT 

C btu 3 

MASS 

EVAPORATED 

Clbm] 


: 6 
s 470 
: 2.44 
: 0.36 
s 28999.99 

1 6.795188 

2 0.0930 
: 0.0636 
t 72 

s • 30000 


2247 
1376.86 
507.1932 
5000 

641.9901 
988.1334 
411.8053 
3.70943E-08 
1000 

12313. 14 
71356.95 
83670.09 


HYDROGEN 

OXYGEN 

988.13 

411.81 

35000.00 

35000.00 

489.00 

489.00 

20. 18 

22. 04 

8. 23 

5.97 

851.S3 

447.96 

0.0285 

0.0226 

2336.53 

891.52 

2. 02 

O. 84 

20.00 

20. OO 

4. 38 

68. 64 

10193.85 

61163.10 

38. 38 

166.41 

188.16 

91.62 

600.00 

600.00 

0. 20 

0. 20 

42. 58 

22. 40 

0. 25 

0. 25 

0.000056 

0.000056 

107.12 

43. 51 

7712.90 

3132.57 

40. 99 

34. 19 
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OXYGEN/HYDROGEN RATIO E. 3 

: 

6 


ISP 

C. 3 

1 

4S6 


ACCELERATION - MAX. 

Eg'a3 

s 

2. 46 


ACCELERATION - MIN. 

Cg'a3 

: 

0.39 


DELTA V - OBTAINABLE 

Cft/aac3 

i 

29000 


HASS RATIO 

C. 3 

: 

6.379754 


STUCTUAL COEFF. 

E. 3 

i 

0. 0988 


PAYLOAD COEFFICIENT 

C. 3 

t 

0.0687 


TIHE, ENDURANCE 

[houra3 

: 

72 


THRUST 

Elbf 3 

: 

30000 


HASSES OF COHPONENTS 

t 




CABIN 

Elba 3 

: 

2247 


LIFE SUPPORT HODULE 

Clba3 

* 

1376. 86 


CARGO HODULE 

Elba 3 

: 

510.5699 


CARGO 

Elba 3 

s 

5000 


SUPPORT STRUCTURE 

Elba 3 

x 

641.3731 


HYDROGEN TANX 

Elba 3 

X 

909.1804 


OXYGEN TANK 

Elba 3 

• 

• 

380.309 


ENGINE 

Elba 3 

t 

1000 


TOTAL BURNOUT HASS 

Elba 3 

X 

12197.84 


HASS OF FUEL 

Elba 3 

X 

65621.38 


TOTAL HASS 

Elba 3 

X 

77819.22 





HYDROGEN 

OXYGEN 

TANK - HASS 

Elba 3 t 

909.18 

380.31 

- HAX. STRESS 

Elbf/in A 23 i 

t 

35000.00 

35000.00 

- HETAL DNSTY 

Elba/ft A 33 


489.00 

489.00 

- PRESS. , BOT. 

£lbf/in A 23 


20.19 

22.13 

- RADIUS 

Eft 3 


8. 01 

5. 81 

- SURFACE AREA 

Eft A 23 


80S.32 

423.62 

- THICKNESS 

Ein3 t 

0.0277 

0.0220 

- VOLUHE 

Eft A 33 * 

2148.97 

819. 87 

- VOLUHE, HTL 

Eft A 33 i 

1.86 

0. 78 

PROP. - VAPOR PRS 

Elbf/ln A 23 


20.00 

20.00 

- DENSITY 

Elba/ft A 33 


4.38 

68. 64 

- HASS 

Elba 3 


9374.48 

56246.90 

- TEHP. 

Erankina3 


38. 38 

166.41 

- HEAT.VAPRZ 

Ebtu/lba3 


188.16 

91. 62 

SURFACE TEHPERATUE 

Cranklne] 


600.00 

600.00 

INSULATION - RHO Elba/ft A 2-in3 i 


0.20 

0. 20 

- HASS 

Elba 3 : 


40. 27 

21. 18 

- THCKNS 

Ein3 : 


0. 25 

0.25 

-K Ebtu-in/hr-ft A 23 i 


0.000056 

0. 000056 

HEAT LEAK - RATE 

Cbtu/hour3 


101.31 

41. 14 

- TOT 

Ebtu3 


7294.41 

2962.40 

HASS EVAPORATED 

Elba 3 


38.77 

32. 33 
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V5TAR Projr^j^ 



10 00SU8 1000 'INITIALIZE VARIABLES 
20 CLSlINPUT*OPTION 1 OR 2 MJOB ... 

30 ON IJOB GOTO 40,210 
. 40 OPEN"O a '» #1, "ROCXETl. OAT* 

30 FOR CONFIG-1 TO 6 * 

60 READ OXH, OXRP, ISP 
70 TEND-72 
' 80 THRUST-30000 

90 TINSULd ) 23 

100 TINSUL(2)-. 23 
110 TINSUL(3)-* 23 
120 PVAPOR(1>-20 
130 PVAPOR(2)-20 
140 PVAPOR(3)-20 
130 GOSUB 2000 
160 GOSUB 3000 'GET MFUEL 

170 PRINT#1, OXH, OXRP, ISP, TEND, THRUST, TINSUL<1), TINSUL < 2 >, TINSUL < 3 >,PVAPOR(1), 
POR < 2),PVAPOR(3), MFUEL 

173 PRINT OXH,OXRP,ISP,TEND,THRUST, TINSUL<1) , TINSUL < 2>,TINSUL(3),PVAPOR <1),PV 

R (2 ) , PVAPOR < 3 > ".FUEL 

160 NEXT CONFIG 

190 CLOSE#1 

200 STOP 

210 GOSUB 7000 ' SORT DATA 
220 OPEN a O a ,#1,*R0CKET3.DAT" 

230 OPEN a I a ,#2, a R0CKET2. DAT a 

240 INPUT#2, OXH,OXRP,ISP,TEND, THRUST, TINSUL <1), TINSUL < 2), TINSUL(3 >, PVAPOR <1) 

POR< 2 > .PVAPOR(3), MFUEL 

230 GOSUB 2000 ' GET FUEL PROPERTIES 

260 GOSUB 3000 * GET PARAMETERS 

270 GOSUB 4000 * PRINT PARAMETERS 

280 IF EOF(2)-0 THEN 240 

290 CLOSE 

300 STOP 


/ 
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iOOO ... INITIALIZE CONSTANTS 

** <3 ? > HVAP0R<3) > HEVAP(3) * MAXSTRS(3) t MTANK(3). PA<3). PTEMP(3>, PVAPOR 
, QDOT(3 ) , QTOT( 3 ) , RADIUS< 3 ) , RFUELO), RINSUL<3>, RMETAL(3>, TANK8<3>, TINSULO), TOl 

1010 DELTAV0-2S0001 
1020 KINSUL<1)>. 000056 
1030 X1NSUL< 2) •.000056 
' 1040 KZNSUL(3)a.0000S6 
1080 HAXSTRSU)■350001 
1090 MAXSTRS<2J>350001 
1100 MAXSTRS(3)*35000i 
1110 HPAY-5000 
1120 PI»3.1415927* 

1130 RMETAL <1)>489 
1140 RMETAL<2)-469 
1150 RMETAL(3)*469 
1160 R1NSUL(1>».2 
1170 RINSUL(2)2 
1180 RXNSUL(3)2 
1190 TANKA(l)••HYDROGEN" 

1200 TANX«<2)••OXYGEN* 

1210 TANK8f3)■*RPJ # 

1220 TOUT(1)*600 
1230 TOUT(2)>600 
1240 TOUT(3)«600 
1250 DIM N(3,3) 

1260 FOR X-l TO 3:F0R Y*1 TO 3:READ N(X, YXNEXT Y:NEXT X 
1270 DATA 8,10, 6 
1280 DATA 10,3,2 
1290 DATA 3,3,2 

1300 DIM THVAPOR (3,6), HVAP0R1 (3,6), PFUEL (3,10), TPFUEL (3,10), RFUEL1 (3, 10), TRFUE 
3,10) 

1310 FOR X*1 TO 3 

1320 FOR Y»1 TO N(X,1):READ PFUEL(X,Y),TPFUEL(X, Y)tNEXT Y 
1330 FOR Y*1 TO N(X,2XREAD RFUEL1 (X, Y), TRFUEL (X, Y) t NEXT Y 
1340 FOR Y-l TO N(X,3)tREAD THVAPOR(X,Y),HVAPQR1(X, Y):NEXT Y 
1350 NEXT X 
1360 RETURN 

1370 DATA 1. 044, 25. 07, 1. 102, 25. 27, 2. 959, 28. 67, 6. 689, 32. 27, 13. 066, 35. 87, 22.'992, 
.47,29.446,41.27,33.903,42.37 

1380 DATA 4. 71, 26. 87,.4. 61, 30. 47, 4. 48, 34. 07, 4. 43, 37. 67, 4. 18, 41. 27, 3. 99, 45. 27, 3. 

, 48. 47, 3. 5, 52. 07, 3. 14, 55. 67, 2. 22, 59. 27 

1390 DATA 32. 27, 898. 83, 35. 87, 860.16, 39. 79, 804. 12, 43. 4, 744. 19, 47, 674. 61, SO. 6, Si 
22 

1400 DATA . 027, 98. 67, . 035, 116. 67, 2. 102, 134. 67, 8. 225, 152. 67, 23. 655, 170. 67, 78. 71 

197. 67, 196. 692, 224. 67, 322. 342, 242. 67, 501. 725, 260. 87, 611. 814, 269. 67 

1410 DATA 85. 65, 97.77,71.2,162.77,.2766,263.67 

1420 DATA 10,91.62,500,91.62 

1430 DATA 10,35,11,35,12,35 

1440 DATA 75,38,75,42,75,46 

1450 DATA 20, 100, 70, 100 

1460 DATA 7,2.21,460 

1470 DATA 7,1000000000000000000,460 

1480 DATA 7,4.25,420 ' 

1490 DATA 6,1000000000000000000,492 
1500 DATA 6,1000000000000000000,486 
1510 DATA 6, lOOOOOOOOOOOOOOOOOO,470 
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2000 »***#*•**#***•**•*•* GET FUEL DENSITIES 
2010 FOR X»1 TP 3 . 

2020 FOR M TO N(X, It / . 

2030 IF PVAFOR(X)<PFUEL(X,Y > THEN 2030 . ' 

2040 NEXT Y 

2030 PTENP(X >-tPFUEL< X, Y-l > ♦ <TPFUEL< X, Y)-TPFUEL < X, Y-l > > * <PVAPOR(X)-PFUEL < X, Y-1) 
/(PFUEL<X, Y)-PFUEL(X,Y-l>> 

2060 FOR Y«1 TO N(X, 2) 

2070 IF PTEMP(X)< TRFUEL<X,Y> THEN 2090 
2060 NEXT Y * 

2090 RFUEL(X)»RFUEL1<X,Y-l)♦<RFUEL1(X, Y)-RFUEL1(X, Y-l))*(PTEMP( X )-TRFUEL< X, Y-l) 
/<TRFUEL(X, Y)-TRFUEL(X,Y-l)) 

2100 FOR Y»1 TO N(X,3) 

2110 IF PTEMP(XX THVAPOR<X,Y> THEN 2130 
2120 NEXT Y 

2130 HVAPOR< X)-HVAP0R1(X,Y-l)♦< HVAP0R1(X, Y)-HVAP0R1(X,Y-l))•(PTEMP < X)-THVAPOR < X 
Y-l)>/<THVAPOR(X, Y)-THVAPOR(X, Y-l)) 

2140 NEXT X 

2130 HVAPOR(1)-HVAPOR(1)/RFUEL(1) 

2160 RETURN 




A 


42 






















3000 '•****#****•*♦♦••*•• FIND CABIN HASS 
3010 HBRNOUT-O 

3020 MCA8IN a 600*lSO*137*1140*100*100 
3030 HBRNOUT-HCABXN 

3040 '•♦***•***•*•*•••*•• FIND LIFE SUPPORT SYSTEM HASS 
3030 HLIFSUP-1073*(33.34/24)»3*TEND 
3060 HBRNGUT a HBRNOUT*MLIFSUP*MPAY 

3070 ... FIND CARGO MODULE HASS 

3080 HCRGM0D a 130*.017*H8RN0UT*AHAX 

3090 HBRNOUT* MBRNOUT+MCRGMOD 

3100 '#*•***♦********••••* fihd tank Basse's 

3110 HFUEL(1)»MFUEL/(1*OXH *OXH/OXRP) 

3120 HFUEL(2)■HFUEL(1)»QXH 
3130 HFUEL(3)*HFUEL(2)/OXRP 
3140 FOR X»i TO 3 
3130 HEVAP«0 

3160 VTANK < X >■<HFUEL< X)♦MEVAP< X))/RFUEL < X) 

3170 RADIUS < X)-<3* VTANK(X)/PI/4)* <1/3) 

3180 AREA(X) a 4*PI»RADIUS(X)*2 

3190 PA(X)*PVAPOR< X > * < RFUEL< X) /12*3)*AHIN* < 2*RADIUS< X)•12) 
3200 TA < X)*PA(X)•RADIUS< X)•12/2/HAXSTRS(X) 

3210 T<X) a TA(X) 

3220 VHETAL<X)"AREA(X)*T(X)/12 
3230 HTANK(X) a VHETAL(X)•RHETAL(X) 

3240 QDOT(X)*KINSUL(X)*AREA(X)* <TQUT<X)-PTEHP(X))/TINSUL(X) 
3230 QTOT(X)«QDOT(X)*TEND 

3260 'PRINT USING"####.####"jflEVAP, HEVAP(X) 

3280 HEVAP < X) a QTOT < X)/HVAPORlX) 

3290 IF A8S(HEVAP-HEVAP(X))>.01 THEN MEVAP a HEVAP<X):GOTO 3160 
3300 HINSUL(X)-RINSUL(X)•AREA< X)»TINSUL(X > 

3310 HBRNQUT*MBRNOUT*HTANK(X)*HINSUL(X)♦HEVAP(X) 

3320 NEXT X 

3330 '•*****••**♦*•*****♦* FIND SUPPORT STRUCTURE MASS 

3340 HSUPSTR•200♦.017*HBRN0UT*AHAX 

3330 HBRNOUT■HBRNOUT♦HSUPSTR 

3360 '•♦*♦•••*••#**♦#*•*•* FIND ENGINE HASS 

3370 HENGINE*1000 

3380 '*•••*♦•••♦**♦*•♦*•♦• FIND BURNOUT HASS 

3390 H BRNOUT a HBRNQUT♦HENGINE 

3400 HTOTAL■HBRNOUT♦HFUEL 

3410 AHAX a THRUST/HBRNOUT 

3420 AHIN-THRUST/HTOTAL 

3430 'PRINT USING”#####.####*;AHAX, AHIN. AC 
3440 IF ABS(AHAX-AC> >.01 THEN AC • AHAX>GOTO 3000 
3430 HRATI0«2. 7182818#*(DELTAVO/32.2/ISP) 

3460 HFUELl a HBRNOUT*(HRATIO-1) 

3470 'PRINT USING"##########.##";HFUEL1{HFUEL;HRATIO 
3480 IF ABS<HFUEL1-HFUEL)>.01 THEN HFUEL*HFUEL1sGOTO 3000 
3490 HRATIO*HTOTAL/HBRNOUT 
3300 DELTAV>32. 2*ISP»L0G<HRATIO) 

3310 EPSILON* < HBRNOUT-HPAY)/(HTOTAL-HPAY) 

3320 LAHBDA«HPAY/(HTOTAL-HPAY) 

3330 RETURN 




























4000 

4010 

4020 

4030 

4040 

4090 

4060 

4070 

4060 

4090 

4100 

4110 

4120 

4130 

4140 

4190 

4160 

4170 

4160 

4190 

4200 

4210 

4220 

4230 


PRXNT#1, 
PRINT#1, 
PRINT#!, 
PRXNT#1, 
PRINT#A, 
PRINT#1, 
PRINT#i, 
PRINT#1, 
PRXNT#1I 
PRINT#!, 
PRINT#!, 
PRINT#!, 
PRINT#1, 
PRINT#1, 
PRINT#1, 
PRINT#!, 
PRINT#1, 
PRINT#1, 
PRINT#!, 
PRINT#!, 
PRINT#!, 
PKINT#1, 
PRINT#!, 
PRINT#!, 


OXYGEN/HYDROGEN RATIO 
OXYGEN/RPJ RATIO 


ACCELERATION - MAX. 
ACCELERATION - MIN. 
DELTA V -OBTAINABLE 
MASS RATIO 
STUCTUAL COEFF. 
PAYLOAD COEFFICIENT 
TIME, ENDURANCE 
THRUST 

HASSES OF COMPONENTS 
CABIN 

LIFE SUPPORT MODULE 
CARGO MODULE 
CARGO 

SUPPORT STRUCTURE 
HYDROGEN TANK 
OXYGEN TANK 
RPJ TANK 
ENGINE 

TOTAL BURNOUT HASS 
HASS OF FUEL 
TOTAL MASS 


C/3 

s'jOXH . 

C. J 

| * 

*•jOXRP 

£. 3 

l"|ISP / 

Cg'a3 

t "t iPRIRT#1 

Lq'mi 

tPRINT#! 

(lt/«#el 

: •;DELTAV 

C. 3 

t • f HRATIO 

t. 3 

1*1:PRINT#1 

C. 3 

: *;PRINT#!, 

£hours3 

: * jTEND 

lib* 3 
• 

* * t THRUST 

(lbs 3 

t»;MCABIN 

Clba3 

t *;HLIFSUP 

ClbaJ 

i* t MCRGMOD 

ClbaJ 

* *>HPAY 

Clba3 ' 

<";MSUPSTR 

ClbaJ 

*”;MTANK(1) 

ClbaJ 

t•jMTANK <2) 

ClbaJ 

1*}MTANK(3> 

ClbaJ 

i *;M£NCINE 

ClbaJ 

t *;MBRNOUT 

Clba3 

i • ;HFUEL 

Clba3 

s•;HTOTAL 
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5000 .PRIHTfi,TAB(2S)»-TAHK # ’ > TAl ’ T \\ PRIMT#i» USING •######♦****♦• 

5010 PRIMT#1,"TANK - MASS ' / 

;-*‘aSKr-“T!‘£iu- STRESS CibN/i«*21,.PRINWl, USING........***- «»' > 

‘,Rliiff‘ 2 ’' "-“n^AL^NSTY cic/lt'3) , • , .PRINT.1, USING-.......».»-.»' > 

■ ’' R " C -*PMSS. . BOT. Iibi/ln-21. ■ , .PRXHT.l, USING-......*•**• *»’ > 

i),PA(2),PA(3) ,**, 

5050 PRINT#1,* - RADIUS £1 ° 

IUSC1),RADIUS(2),RADIUS(3) £1^*23 

5060 PRIMT#lp * * SURFACE AREA tit 2 

A(1># AREA<2)» AREA<3) £xa j 

5070 PRINT#1," “ THICKNESS 

<l>#T<2),Tt3> tft R 3J 

5060 PRINT#i,“ ' VOLUME 3 

NK<1>,VTANK(2>,VTANK<3i 
5090 PRINT#!,• * VOLUME, MTL tit 

TAL(l), VMETAH2), VMETALt3^ OBfl ci b x/in*21 * * t :PRIHT#1, UPING*##*#*#**##* **' f * 


,•,:PRINT#1, USING•##########*#*"* 
, *, jPRINT#1, USING*#####***##. **' '* 
x• ; x PRINT#1, USING*#######** *#**" 

*• { tPRINT#i,USING*######*#**-**’ ; 
i • i -.PRINT#!, USING*##»*******♦ **’ ; 


£15«/*t*3J s * t sPRINT#l, USING*####****#*•**" ; 1 
Llbmi s"j sPRINT*1, USING*##*#*#*###•#*"' 1 

crmnkin#! s * ,* sPRINT#!, USING*#####*##**- **' ;f 



; r 


5100 PRINT*!,* PROP.- VAPOR PR» 

POR11),PVAPOR(2),PV APOA< 3 > 

5110 PRINT#1,* “ DENSITY 

EL(1),RFUEL < 2 >,RFUEL < 3) 

5120 PRINT#!,* “ M* 53 

EL( 1), MFUELv2>» MFUEH3) 

5130 PRINT#1,* • TEMP. 

"““i^SAT. VAPRZ Cb.u/ib.l •’>-PRINT.I. USING .. 

SISO^RINT.!^■SURF»CE°TE»PNRAVUE Cr.nNAn.i .PRINT.l.USING-...‘^ 1 

Jllo' PRINTill * INSULAYiUN - HNO Cltom/*«~2-in J s * f sPNxNT#!, UsINO 

SUH 1 / , RINSUi»v 2 j , RINSULvwK « :PRINV#1, USING'#****###*** **" **■ 

5170 PRINT#1. * _* nA * a LJ - DB,J 

SUL<1>, M1NSUH2), MINSULv , • . ,PRINT#!, USING*####*##*#** **" • 1 

5160 PRINT#!,* THORNS ixnj 

SUL(1>,TINSUL(2),TINSULvS> .„ tu . 1B/hr . <t *2]:* i JPkINT#l, USING*######.####»#*;" 
5190 PRINT#1,* * £totu in/nr 

"* TE tbtu/hourl t"(iPRINT.i» USING-........ 

T(1>» ODOT<2>,QDOT<3> ,•,»PRINT#1,USING*####***###**** 

5210 PRINT#!, * * TOT 

T11>, OTOT12), UTOT13 > D c Ato- j ,-;sPRINT#!, USING * ##»###*##*.#*"J« 

5220 PRINX*i» * EVAPURATED ixobm 

AP(1),MEVAP12>,MEVAP<3) 
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6000 H ANV6<1)■•INJECT HOHHAN TRANSFER/" 

6010 HANV*<2>*"CIRCULARIZE 

6020 ffANV«<3)-”D£-CIRCULARIZE >'• 

6030 HANV*<4>-'CIRCULARIZE " 

6040 DELTAV<l>»d929. 87 
6090 DELTAV(2)*9741.19 
6060 DELTAV<3>-9741.19 
r 6070 DELTAV(4)*6929. 67 
6060 T-HFUEL-ISP/THRUST/60 
6090 PRINT#1, "TINE OF BURN CminutM] r";T 
6100 HIN1T-MT0TAL 
6110 PRINT}FOR X»1 TO 4 

6120 MFINAL-HINIT-2. 72*(-DELTAV(X)/32. 2/ISP) 
6130 HFUEL-HINIT-HFINAL 
6140 T-NFUEL-ISP/THRUST/60 

6190 PRINTFl," - ■ j HANV*<X)j"i";T, DELTAV(X > 

6160 HINIT-HINI7-HFUEL 
6170 NEXT X 

6160 FOR X»1 TO 19:PRINT#1,•-:NEXT X 
6190 RETURN 
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4 

* * , * * * ~ 

7°i° DIH A<H # <H^o CKBT1 DaT . ^ 

SS-b-: .2 .-oc“ ts - dat v 

%£ r-iir 0 A(H( 

^rp p s«r 0 v 050 

7090 IF 0 #**’ > A THCM A ‘ U<X>, ' , ' X £(B)1 F1-F<M>'01- 6< " ) ' H1 “ B< " > '“' I ‘"’ ! 

4t?S t 5 > • ? i:»‘ ( Si. cx- c < « » . OX-0 <«*..“■ 

7i30 A(Y>»Ai* B(Y> Bi 

7150 PRIHT#2,A<«>»»' n, * w 

71^0 them 70d0 

7170 IF * <-M THEN 
7ISO CLOSE 1 
7190 RETURN 
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« • • 

MOO OPEN "OS #1# "FUEL*DAT* / 

8010 FOR X>1 TO 3 . ' 

8020 FOR Y«1 TO N(X, IMPRINT#!, PFUEUX, Y), TPFUELtX, YMNEXT Y 
8030 FOR Y*1 TO N < X, 2 M PRINT#!, RFUEL1< X, Y >,TRFUEL< X, Y MNEXT Y 
8040 FOR Y»i TO N ( X, 3) * PRINT#!, THVAPOR < X, Y), HVAPORIX, Y M NEXT Y 
8030 NEXT X 
8060 CLOSE#! 




















Temperature at boiling point vs. Pressure 
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Pressure (PSI) 
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Density of liquid at boiling point vs. Temperature 
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Temperature (Rankine) 
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Heat of vaporization vs. Temperature 
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